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Doping of a two-orbital chain with mobile S = 1/2 Fermions and strong Hund's rule couplings 
stabilizing the S = 1 spins strongly depends on the presence of a level difference among these 
orbitals. By DMRG methods we find a finite spin gap upon doping and dominant pairing 
correlations without level- difference, whereas the presence of a level difference leads to dominant 
charge density wave (CDW) correlations with gapless spin-excitations. The string correlation 
function also shows qualitative differences between the two models. 
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Hole doping of strongly coupled spin-liquid ground- 
states is of great importance for the understanding of 
the Mott-Hubbard transition. A very special spin-liquid 
is the antiferromagnetic (AF) S = 1 Heisenberg (HB) 
chain, which shows a finite spin gap 1 ) and a hidden topo- 
logical order, the string-order parameter. 2 ) Only recently 
it has become possible to study experimentally the ef- 
fect of doping in a spin-1 chain with mobile holes in 
the system Y2- x Ca x BaNi05. 3,4 ) In this material, the 
5 = 1 spins are formed by a strong ferromagnetic (FM) 
Hund's coupling Jh between the two active orbitals of 
Ni 2+ . Theoretically, a number of interesting questions 
arises by doping of a S = 1 chain with mobile holes. 
The competition between FM order induced by the dou- 
ble exchange mechanism 5 ) and AF order can result in 
completely different magnetic properties, the spin gap 
can immediately be destroyed upon doping or might per- 
sist, and finally which correlation function characterized 
by the single correlation exponent K p dominates in the 
thermodynamic limit. 

In a previous paper, 6 ) we have investigated the prop- 
erties of a doped spin-1 chain with a level difference be- 
tween the two electrons forming the S — 1 spin. Such 
a situation is likely in a system where the S = 1 spins 
are formed by two electrons in different orbitals of the 
same ion, as it is the case in Y 2 _ x Ca x BaNi0 5 . 7 - 9 ) Here 
we will concentrate on the effect of a level difference, 
by comparing the above case with a system with elec- 
trons in two equivalent levels. The latter situation is 
more likely in a ladder-system. A few years ago, this 
problem has been studied by Fujimoto and Kawakami 
in a weak coupling approach, where two ferromagneti- 
cally coupled Hubbard-chains have been investigated and 
compared to a Hubbard chain ferromagnetically coupled 
to a Heisenberg-chain. 10 ) In their analysis it is found 
that the model with two equivalent chains remains gap- 
ful, with a finite string-order parameter and dominant 
CDW correlations (K p — 1/2). In contrast, the model 
with a level-difference becomes gapless upon doping and 



shows two gapless spin-modes. Similar conclusions have 
been drawn by Nagaosa and Oshikawa in a semi-classical 
approach with two coupled Hubbard chains, 11 -* and it is 
pointed out that for low hole doping, the system remains 
gapful for equivalent doping on the chains irrespective of 
FM or AF couplings on the rungs, whereas the system 
becomes gapless upon doping for different hole concen- 
trations on both chains. However, since both of these 
studies are based on weak-coupling or mean-field ap- 
proaches, it is not clear whether the results are valid 
in this one-dimensional, strongly correlated system and 
direct numerical evidence is needed. In this Letter we 
would like to investigate and compare doping of a spin- 
1 chain with and without level difference in an unper- 
turbative manner free of approximations by the density 
matrix renormalization group method (DMRG). 12) This 
method is ideal for the study of groundstate energies and 
equal-time correlation functions, and we determine the 
correlation exponent K p and the string order parame- 
ter. Additionally we also calculate the magnetic sus- 
ceptibility at finite temperatures by the thermal DMRG 
(TDMRG) method. 13 ' 14 ) 

The first model we consider is a model with equal par- 
ticle concentrations on both legs, and we call this the 
symmetric model in the following. The S = 1 spins are 
formed by two ferromagnetically coupled t-J chains with 
additional AF couplings between next-nearest neighbor 
sites on opposite legs. The individual terms which define 
the model are given by hopping along the chain i 
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where the projection operator V prohibits doubly occu- 
pied sites, and ct . a is the particle creation operator on 
rung j and leg i with spin a. Further the AF couplings 
J > on the chain i between nearest neighbors and the 
diagonal couplings Jd > are given by 
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3,1 V ' 

where the notation is standard, the indices are the same 
as above, and 5jj denotes Kronecker's delta function. 
The strong Hund's coupling Jh < acts on the rungs 

i?FM = Jh ^ Sj t iSj t 2- 

j 

Finally our Hamiltonian reads 

H sym = H FM + tfdiag + J2 ( H * + ^kin) • 

i 

The second model contains a level difference between the 
two electrons forming the S = 1 spins and the lower band 
is localized. We will refer to this model as the asym- 
metric model in the rest, and it consists of a t-J chain 
ferromagnetically coupled to a S = 1/2 HB chain, with 
additional AF couplings between next-nearest neighbors 
on opposite legs. The Hamiltonian is defined as 

-ffsym = HFM + Hding + + H^. 

For both models we restrict ourselves to parameter values 
with J = J d and \J H \ » t, J, J d , and if not otherwise 
mentioned we set —Jh = Wt = 20J = 20Jd- In the 
DMRG calculations we use systems sizes of up to 2 x 256 
sites and up to 1300 states per system- and environment- 
block. In the TDMRG we keep 80 states per block and 
use finite Trotter-time steps of At = 0.2t. 

We have reported on the asymmetric model in a previ- 
ous paper, 6 - ) and we briefly summarize the main results 
here. In that system, each mobile hole creates a small 
FM cloud by the double exchange mechanism (polaron) . 
There is a very weak AF interaction among the polarons 
giving the lowest lying, gapless spin-excitations. How- 
ever, the polarons are only a weak perturbation on top 
of the underlying spin-liquid of the Haldanc chain, and 
there is a second energy scale in the spin-sector of the 
order of the Haldane gap of the undoped system. This 
energy scale shows up in exponentially decaying spin- 
correlations at short to intermediate distances and a fi- 
nite string-order. The dominant correlation in the ther- 
modynamic limit is given by 2/cp and Akp CDW order 
and the correlation exponent is K„ w 0.51. 

In the following we will perform a similar analysis for 
the symmetric model and compare these findings with 
the asymmetric model. We start with the investigation 
of the spin gap. At half filling and Jh ^> J,Jd, both 
models can be mapped to the Haldane S = 1 chain 
with effective couplings J = J^g m in the symmetric case 
and JJg = 3/4 J in the asymmetric case. By finite size 
scaling we have numerically determined the spin gap 
from A s = limp^o A S (L; N = Ln), where A S (L;N — 
Ln) = A S {L-N) = E {L-N;S Z = 1) - E (L;N;S Z = 
0) and Eq(L; N; S z ) is the groundstate energy of the 
system with N particles on L sites and the total spin 
component along the z-direction is S z . The results are 
in excellent agreement with the above mapping and we 
find A^™ w 0.41(l)J c s y. m = 0.205(5)i for the symmetric 
model and A^ in w 0.41(1) J c a | w 0.154(5)i for the asym- 
metric model. In agreement with refs., 10 ' we find com- 
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Fig. 1. Temperature dependencies of the magnetic susceptibility 
\ at different particle densities. Symmetric model (a) with J = 
Jd = t/2 and Jjj — > oo, and asymmetric model (b) with —Jh = 
10* = 20 J = 20 J d . 



pletely different behavior for both models upon doping, 
as can be seen in Table I. While the symmetric model 
retains a rather large spin gap which decreases slightly 
upon doping (n is the particle density per single site), 
it is destroyed already for the smallest hole density of 
n c h = 0.0625 holes per site of the conduction band for 
the asymmetric model. 

At finite temperatures near the gap value T oc J/2, the 
magnetic susceptibility x of both systems is strongly en- 
hanced upon doping, and we show the TDMRG results 
in Fig. 1, where we have set Jp — > oo for the symmetric 
model for numerical reasons. This strong enhancement 
of x is caused by the creation of a small FM cloud around 
each mobile hole by the double exchange mechanism, as 
the holes can gain kinetic energy by a FM alignment of 
the neighboring spins. At still lower temperatures, the 
magnetic susceptibility is suppressed for the symmetric 
model due to the spin gap and we find x ~~ * for T — > 0. 
In contrast, \ is finite and rather large at T = for 
the asymmetric model, indicating the formation of larger 
FM moments due to the double exchange mechanism (see 
ref. 6). However, the groundstate is spin-singlet for both 

Table I. Finite size scaling of the spin gap A sp i n for — Jh = 10t = 
20 J = 20 Jd and different values of the hole density n/, = 1 — n for 
the symmetric model and hole density in the conduction band n c h 
for the asymmetric model. Also listed is the spin-spin correlation 
length £. 
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0.205(5)t 


6.05(5) 


0.1504(5)* 


6.05(5) 


0.0625 


0.148(5)t 
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0.125 


0.120(5)* 


6.7(2) 





11.2(2) 
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Fig. 2. Spin-spin correlations (a) for the symmetric and asym- 
metric model with — Jjj = lOt = 20 J = 20 Jd, and string-order 
g(x) for both models with same parameters (b). 



Fig. 3. Charge density for the symmetric model (a) at n = 7/8 
and the symmetric model (b) at n c = 7/8, the parameters for 
both models are — Jh = lOt = 20 J = 20 and the length 
is L = 128. In the insets we show the corresponding Fourier 
transforms n(k). 



models, as we have tested with the DMRG method by 
the calculation of a(i,j) = 5?5|- 1/2(5/" ST") which van- 
ishes for a rotationally invariant groundstatc and the re- 
sult is zero within the numerical precision of the DMRG 
for this quantity for both models (\a(i, j)\ < 10~ 6 Vi, j). 

Complementary information about the magnetic prop- 
erties is provided by the spin-spin correlations shown in 
Fig. 2(a). At half filling we find exponentially decay- 
ing correlations as expected for gapful systems and the 
correlation length is £ s=s 6.03. 15 - ) Also at finite dop- 
ing, the spin-spin correlations seem to decay exponen- 
tially for both systems. We have determined £ by fit- 
ting to SfSf +x oc cos(2kpx)e~ x ^ and the results are 
listed in Table I. In the case of the symmetric model, 
the correlation length increases only slightly upon dop- 
ing to £ sym « 6.7 ± 0.2 for n/, = 1/8, while the increase 
to £ as w 11.2 is considerably larger for the asymmetric 
model at the same density. In fact we expect a crossover 
to power-law decay at large distances in the latter case 
(see ref. 6). More insight about the spin-configuration 
can be gained from the string-order parameter g(x) = 

((Ei=l,2^«,i) (rife=xo+l,J=l,2 e l7rSfc ' 3 ^ (El=l,2^,l))> 

which reveals the hidden Z-i x symmetry of the Hal- 
dane 5=1 chain and quickly approaches the value of 
g(x) w —0.374 for x » 1 in the undoped case. 2 ) From 
Fig. 2(b) it can be seen that \g{x)\ is reduced upon dop- 
ing but remains finite for both models. A finite string 
order parameter for the symmetric model has also been 
obtained in ref. 10. Compared to the symmetric model, 
the reduction of \g(x) | upon doping is roughly three times 
larger for the asymmetric model. In addition, the oscil- 
lations in |<7(x)| are a strong amplification of the cor- 
responding Friedel oscillations in the charge density in- 
duced by the open boundary conditions. In fact, for the 
asymmetric model, we again expect a crossover of \g(x)\ 



to power-law decay at large distances after a sufficient 
number of oscillations. By comparison, \g(x)\ is con- 
stant within 1.2% for the symmetric model apart from 
boundary effects, and we conclude that its spin structure 
closely resembles that of the undoped case. 

Let us turn to the charge sector next and start with the 
above mentioned Friedel oscillations in the charge den- 
sity. In their Fourier transforms we only find one peak at 
k = 2nir for the symmetric model. This corresponds to 
the 2kp CDW fluctuations, since the lower band is fully 
occupied and hence there are 2(1 — n) holes in the upper 
band, giving 2kp = (1 + 1 — 2(1 — n))-n = 2nn for a large 
Fermi volume. In the case of the asymmetric model we 
also find only one peak at k — n^n for the smallest hole 
densities n c h = 1/16, and two peaks at k = n c y jx and k — 
2n^7r for the larger hole doping n c h = 1/8. This is con- 
sistent with 2k p and Akp CDW fluctuations for a large 
Fermi volume including the electrons of the lower orbital 
with 2kp — (n c +l)7r in that case. Having identified these 
2kp and Akp fluctuations, we can determine the correla- 
tion exponent K p by fitting to the Friedel oscillations of 
an impurity potential. In the gapful case they are given 
by 6 n (x) oc C\ cos{2k fx)x~ k i> + C2 cos(4fci?x)x _2 - ft - p and 
S n (x) cx C\ cos(2k F x)x- ( - 1+K »V 2 + C 2 cos(4fc F a;)a;- 2 ^ 
in the gapless case. 16 ) The correlation exponents thus 
obtained are K p w 1.5 ± 0.05 for the symmetric case at 
rih = 7/8 indicating dominant superconducting pairing 
correlations, and K p w 0.51 ± 0.05 for to asymmetric 
model at n c h = 1/8, giving dominant CDW correlations. 
For both cases, the sites near the boundaries need to be 
discarded for the fit because of trapped states, and the 
uncertainty stems from the fitting ambiguity. In the sym- 
metric model hole pairs are formed on the rungs in order 
to gain the strong Hund's rule coupling Jf- The pair 
binding energy A pa i r w 2.29i in the low-doping region 
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Fig. 4. Singlet pairing correlations Pi(f)P^ +x (f) with Pj{f) = 

^75( C !,T C !+/,| _ c [| c !+/,t) for ( a ) s y mmetric model at n = 7/8 
and (b) asymmetric model at n c = 7/8 and —Jh = 104 = 20J = 
20J<j for both models. 



obtained from A pa ; r = 2E\ — Eq — Ei is correspondingly 
large, where E n is the groundstate energy with n holes. 
This pair-binding mechanism does of course not work for 
the asymmetric model, but we still find a positive pair- 
binding energy A pa j r « 0.016i in that case. Pair forma- 
tion in that context has also been found in ref. 9. From 
a second structure that develops in the hole-pockets at 
larger doping, we expect that these pairs get spatially 
more extended upon doping. 

The pairing correlations defined as Pj(f )P- +X (i) with 

P i(?) = 75( c i,T c !+f.l T c [| c i+f.T) allow an independent 
calculation of the correlation exponent K p , and provide 
further information on the form- factor f of the pairs. 
From the numerous possibilities of the form-factor f for 
the symmetric model, we restrict the discussion to those 
with the largest amplitudes in the following. At very 
short distances d < 5 triplet pairing correlations for 
rung pairs with f = (0, 1) have the largest amplitudes, 
but as the triplet correlation functions decay exponen- 
tially, the amplitudes for singlet pairs formed on the 
rungs with an extension of the pair over two to three 
rungs become larger at distances d > 5. We show the 
results for the singlet pairs with the largest amplitudes 
in Fig. 4. These correlations decay as Pi{f)P} +x {f) oc 
x ~ x / k p for a gapful Luther-Emery liquid, 17 -* and surpris- 
ingly there are almost no 2k f fluctuations for the largest 
pairing correlations despite the Friedel oscillations in the 
charge density. The correlation exponent determined 
from the fit K p w 1.55 ± 0.05 is in excellent agreement 
with the result obtained from the charge-density. For 
comparison, the correlation exponent obtained by weak- 
coupling is K p = 0.5. 10 ) The form factor of the pairs 
with largest amplitudes f = (1,1), (2,1) and (3,1) are 
consistent with a 4- 9 symmetry analogue for a ladder 



with vanishing f = (0, 1) and (1,0) amplitudes. For the 
asymmetric model we obtain K p w 0.51 ± 0.05 for a 
gapless Tomonaga-Luttinger liquid with Pi{f)P} +x {f) oc 
A Q ln(x)- 1 - 5 x- 1 - 1 / K "+A 2 cos(2fc F x)x-^- 1 /^pi8) by si _ 
multaneously fitting to / = 2, 4, 6, and 8, again in excel- 
lent agreement with the previous fit. 

In conclusion we have studied the effect of a level dif- 
ference on doping of two-orbital chains. In addition to 
a finite spin gap for the symmetric model and gapless 
excitations for the asymmetric model also obtained by 
weak-coupling theories, 10, we find a new difference in 
the dominant correlation functions, which consist of pair- 
ing correlations with K p w 1.5 for the symmetric model 
and CDW correlations with K p w 0.5 for the asymmetric 
model. Further the string-order parameter remains finite 
for the symmetric model, however strongly reduced by 
the polarons in the asymmetric model. 
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